SPECIFICATION 



TO ALL WHOM IT MAY CONCERN: 

BE IT KNOWN THAT WE, Kosei Takiishi, a citizen 
of Japan residing at Yokosuka-shi , Kanagawa Japan and Tokio 
Taga, a citizen of Japan residing at Yokosuka-shi, Kanagawa 
Japan have invented certain new and useful improvements in 

A MULTI-PATH GENERATING APPARATUS, A MULTI-PATH FADING 
SIMULATOR, AND A MULTI-PATH GENERATING METHOD 



Of which the following is a specification:- 
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TITLE OF THE INVENTION 

A MULTI-PATH GENERATING APPARATUS, A 
MULTI-PATH FADING SIMULATOR, AND A MULTI-PATH 
GENERATING METHOD 

5 

BACKGROUND OF THE INVENTION 

1. Field of the Invention 

The present invention generally relates to 
a time-varying multi-path generating apparatus, a 

10 multi-path fading simulator employing the time- 
varying multi-path generating apparatus, and a 
multi-path generating method for simulating multi- 
path fluctuations in radio communications systems, 
such as mobile communications systems; and 

15 specifically relates to a time-varying multi-path 

generating apparatus, a multi-path fading simulator 
employing the time-varying multi-path generating 
apparatus, and a multi-path generating method 
wherein time-varying multi-paths are generated by 

20 adding time-varying propagation path fluctuations to 
parameters such as arrival direction characteristics, 
phase characteristics, time delay characteristics, 
and power level characteristics of each propagation 
path corresponding to input parameters. 

25 2. Description of the Related Art 

Fig. 1 is a conceptual diagram showing the 
case wherein a radio propagation path is provided 
from a base station antenna to a mobile station, the 
radio propagation path often reaching the mobile 

30 station through multiple propagation paths 

established by reflection and diffraction due to an 
obstacle such as a building. On the right-hand side 
of Fig. 1, an example of a time delay profile is 
shown. In the time delay profile where five 

35 propagation paths are shown, the horizontal axis 
represents the propagation delay time (hereafter 
called time delay) of the radio propagation paths 
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(hereafter called propagation paths) that arrive at 
the mobile station, and the vertical axis represents 
received power. The time delay and the received 
power may be expressed in either of absolute values 
5 and relative values. 

The propagation paths (1), (2), (3), (4), 
and (5) shown in the delay profile are regarded as 
plane waves. A propagation path that has the 
smallest time delay is the propagation path (1) that 

10 reaches the mobile station via the shortest path 
from the base station, and the propagation paths 
with greater time delays than the propagation path 
(1) reach the mobile station after reflection and 
diffraction by buildings, mountains, etc. As for the 

15 propagation paths with the greater time delay, since 
the propagation loss increases as the propagation 
distance becomes greater with the number of times of 
reflection and diffraction increasing, received 
power tends to become smaller. 

20 Fig. 2 is a graph that shows an example of 

the delay profile measured with a mobile station 
moving along a road (non-patent reference 1) . It is 
clearly shown that the time delay and received power 
of a propagation path are changing with the movement 

25 of the mobile station. 

When a mobile communications system is 
designed, a simulator for simulating the 
transmission quality of the space between the base 
station and the mobile station is required in order 

30 to determine whether the system can allow certain 
propagation path fluctuations, and in order to 
determine optimal system parameter values. A fading 
simulator is often used for evaluating the 
transmission quality of mobile communications 

35 systems. The fading simulator applies fading to a 

modulated signal that is provided by a transmitter, 
and outputs the faded modulated signal to a receiver. 



In that case, the fading simulator superimposes a 
modulating signal on a carrier propagation path that 
has a time delay, the time delay being specified by 
a simulator operator according to a profile model 
such as the Vehicular-B model (non-patent reference 
2) recommended by ITU (International 
Telecommunication Union), etc., and applies 
Nakagami-Rice fading or Rayleigh fading; or, 
alternatively, the delay profile is calculated by 
ray tracing, wherein modulating signals are 
superimposed on carriers having different time 
delays, and fading that is. the same as above is 
applied. 

In conventional mobile communications 
systems with a relatively low transmission speed, it 
is sufficient to consider receiving power level 
characteristics (e.g., Rayleigh fading model; non- 
patent reference 3), and delay characteristics (e.g. 
equivalent transmission-line model; non-patent 
reference 4). For this reason, it is not necessary 
to generate propagation characteristics of a 
propagation path that properly simulates the real 
environment. Accordingly, fading simulators are 
designed to provide power level fluctuations that 
follow the Nakagami-Rice distribution or Rayleigh 
distribution that defines propagation path 
fluctuations as a vector sum of two or more 
propagation paths without taking correlations 
between propagation paths into consideration. 

[Non-patent reference 1] 

Masaaki Shinji, "Radio propagation path 
propagation in wireless communications", 1992, page 
210 

[Non-patent reference 2] 

ITU-R M.1225, "Guidelines for evaluation 
of radio transmission technologies for IMT-2000", 
1997 
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[Non-patent reference 3] 

Jakes, W.C.Jr. and ed. : "Micropropagation 
Path Mobile Communications", John Wiley & Sons, Inc., 
New York, 1974 
5 [Non-patent reference 4] 

Y. Karasawa, T.Kuroda, and H.Iwai: "The 
equivalent transmission-path model, - A Tool for 
Analyzing error floor characteristics due to inter- 
symbol interference in Nakagami-Rice fading 
10 environments", IEEE Trans. Veh. Technol . 46 [1] pp. 
194-202, 1997 

[Non-patent reference 5] 

Yamada, Tomisato, Matsumoto: "The system 
evaluation of a spatial-temporal equalizer using 
15 real propagation data" Shin-Gaku Giho A-P 2000-96 
(2000-10) 

[Non-patent reference 6] 

S. Ichitsubo, K. Tsunekawa, and Y. Ebine: 
"Multipath Propagation Model of Spat io-Temporal 

20 Dispersion Observed at Base Station in Urban Areas", 
IEEE Journal on selected areas in communications, 
Vol.20, No. 6, August 2002'. 

[Problem(s) to be solved by the Invention] 
Indispensable technology for the 4G 

25 (fourth generation) mobile communications 
(broadband- transmission) includes adaptive 
equalization technology for removing an unnecessary 
propagation path on the delay time axis, and 
adaptive array antenna (AAA) technology for 

30 separating an unnecessary propagation path on the 
space axis (the radio propagation path arrival 
direction) . Research and development towards the 
realization of these technologies are being 
energetically furthered (non-patent reference 5) . 

35 Further, where the transmission speed (bit rate) is 
high, an incoming propagation path can be separated 
into multiple propagation-delay propagation paths (a 
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single propagation path, or multiple propagation 
paths compounded) . For this reason, the influence of 
fading that is a problem in conventional 
transmission systems at a low transmission speed can 
5 be mitigated in high-speed transmission, and 

communication quality can be improved. Accordingly, 
in order to perform highly precise evaluation of 
transmission characteristics, it becomes important 
that propagation characteristics be individually 

10 reproducible for each of such multiple propagation 
paths between the base station and the mobile 
receiving station. 

The propagation characteristics of a 
propagation path to be considered include 

15 propagation-delay characteristics and arrival 
direction characteristics, and analysis of the 
characteristics and modeling thereof are being 
studied (non-patent reference 6) . For example, it is 
known that an exponential-function model reasonably 

20 approximates the propagation-delay characteristics, 
namely, the received power exponentially decreases 
according to the time delay (non-patent reference 3). 
However, in model approaches such as described above, 
the time variation characteristics of the 

25 propagation path are defined as a mere 

superimposition of either the Nakagami-Rice fading 
or the Rayleigh fading fluctuations, where 
fluctuation characteristics of the multiple 
propagation paths are simply added. 

30 Furthermore, in order to evaluate a 

communications system using adaptive equalization 
technology and AAA technology, a propagation path 
model needs to consider all of the arrival direction 
characteristics, delay time characteristics, and 

35 receiving power level characteristics. 

As a factor of the power level fluctuation 
of each propagation path, shadowing (blocking) of 
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the propagation path by the presence and movement of 
an obstacle such as a car and a human body is to be 
considered. In addition, when the mobile station 
moves, a building and the like may interfere with 
5 the propagation path. For this reason, evaluation of 
each propagation path has to take the power level 
fluctuation caused by these obstacles into 
consideration; however, conventional fading 
simulators do not consider these matters. 

10 Further, as for the propagation paths 

received by the mobile station in directions that 
are approximately the same, there must be meaningful 
correlation between the propagation paths, since the 
same obstacle would interfere with each of the 

15 propagation paths. However, a fading simulator that 
considers the correlation characteristics of such 
propagation paths is not yet offered. 

SUMMARY OF THE INVENTION 

20 Accordingly, it is a general object of the 

present invention to provide a time-varying multi- 
path generating apparatus, a multi-path fading 
simulator using the time-varying multi-path 
generating apparatus, and a multi-path generating 

25 method that substantially obviate one or more of the 
problems caused by the limitations and disadvantages 
of the related art. 

Features and advantages of the present 
invention are set forth in the description that 

30 follows, and in part will become apparent from the 

description and the accompanying drawings, or may be 
learned by practice of the invention according to 
the teachings provided in the description. Objects 
as well as other features and advantages of the 

35 present invention will be realized and attained by a 
time-varying multi-path generating apparatus, a 
multi-path fading simulator using the time-varying 
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multi-path generating apparatus, and a multi-path 
generating method particularly pointed out in the 
specification in such full, clear, concise, and 
exact terms as to enable a person having ordinary 
5 skill in the art to practice the invention. 

To achieve these and other advantages and 
in accordance with the purpose of the invention, as 
embodied and broadly described herein, the invention 
provides a time-varying multi-path generating 

10 apparatus for generating time-varying propagation 
paths, each of the time-varying propagation paths 
being generated based on the time variation 
characteristics (time -varying mult i -path 
fluctuation) of propagation path parameters such as 

15 power level characteristics, phase characteristics, 
time delay characteristics, and arrival direction 
characteristics, corresponding to simulation 
conditions that are input with the correlation 
characteristics of the propagation paths being taken 

20 into consideration; a multi-path fading simulator 
using the time-varying multi-path generating 
apparatus, and a multi-path generating method. 
[Means for Solving the Problem] 
An aspect of the present invention 

25 provides a time-varying multi-path generating 

apparatus that includes a parameter control unit, a 
data storage unit, a random number generating unit, 
and a propagation path generating unit. 

Another aspect of the present invention 

30 provides a time-varying multi-path generating 

apparatus that includes a parameter control unit, a 
data storage unit, a random number generating unit, 
a time-varying function generating unit, and a 
propagation path generating unit. 

35 Another aspect of the present invention 

provides a multi-path fading simulator that includes 
the time-varying multi-path generating apparatus, a 
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propagation path output unit, and a signal 
synthesizing unit. 

Another aspect of the present invention 
provides a time-varying multi-path generating method 
5 that includes steps such as a step for storing 

various propagation path generation parameters and 
data files, and so on, such that multiple 
propagation paths are generated, wherein the steps 
include a step for the parameter control unit 

10 reading a file prepared for generating constant 
parameters of the time-varying functions. 

Another aspect of the present invention 
provides a time-varying multi-path generating method 
that includes steps such as a step for storing 

15 various propagation path generation parameters and 
data files, and so on, such that multiple 
propagation paths are generated, wherein the steps 
includes a step for the parameter control unit 
specifying a calculation model based on a 

20 propagation path generation parameter file. 

[ Function] 

Accordingly, the present invention does 
not simulate propagation path fluctuation based on 
Nakagami-Rice fading, Rayleigh fading, etc., that 
25 are' fluctuation characteristics of a synthetic 

result of each propagation path. But rather, the 
present invention considers propagation path 

fluctuation generated by shadows due to obstacles v 
including moving objects such as vehicles and human 

30 bodies, and buildings that are actually present in 
real propagation paths. Further, the present 
invention also considers the correlation between 
spatial propagation paths. In this manner, the 
present invention provides the time-varying multi- 

35 path generating apparatus, the multi-path fading 
simulator using the time-varying multi-path 
generating apparatus, and the multi-path generating 



method for improved simulation of multi-path fading. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 shows conceptual examples of 
5 multiple propagation paths of a radio propagation 
path provided from a base station to a mobile 
station through reflection and diffraction, and a 
time delay profile of radio propagation through each 
propagation path; 
10 Fig. 2 is a graph showing an example of 

measurement of the time delay profile when the 
mobile station is in motion; 

Fig. 3 is a block diagram showing a 
configuration example of a time-varying multi-path 
15 generating apparatus according to an embodiment of' 
the present invention; 

Fig. 4 is a block diagram showing a 
configuration example of a multi-path fading 
simulator according to the embodiment of the present 
20 invention; 

Fig. 5 is a 3-D graph showing a phase 
fluctuation when the mobile station is in motion 
according to the embodiment of the present 
invention; 

25 Fig. 6 is a flowchart showing operations 

of the time-varying multi-path generating apparatus 
according to the embodiment of the present 
invention; 

Fig. 7 shows graphic examples of initial 
30 values of the propagation path parameters generated 
according to the embodiment of the present 
invention; 

Fig. 8 shows graphic examples of a time- 
varying function according to the embodiment of the 
35 present invention; 

Fig. 9 shows another set of graphic 
examples of the time-varying function according to 
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the embodiment of the present invention; 

Fig. 10 shows graphs for explaining 
propagation path fluctuation output according to the 
embodiment of the present invention when shadow 
5 timing is overlapped; 

Fig. 11 is a flowchart showing operations 
of a time-varying function generating unit according 
to the embodiment of the present invention; 

Fig. 12 is a graph of an example for 
10 calculating a power level fluctuation due to shadow; 

Fig. 13 is a graph showing comparison 
results of the power level fluctuation calculation 
with measured values concerning the shadow; 

Fig. 14 shows graphic examples of 
15 receiving power levels when the power level 

fluctuation due to shadow is provided to each 
propagation path; 

Fig. 15 is a block diagram showing a 
configuration example of a propagation path output 
20 unit according to the embodiment of the present 
invention; 

Fig. 16 is a block diagram showing a 
configuration example of a signal synthesizing unit 
according to the embodiment of the present 
25 invention; 

Fig. 17 is a graph that shows an example 
of the results that the fading simulator outputs 
according to the embodiment of the present invention 
with the Rayleigh distribution being added; and 
30 Fig. 18 is a graph that shows an example 

of the measurement result of the synthesized 
receiving level with the Rayleigh distribution being 
added. 

35 DESCRIPTION OF THE PREFERRED EMBODIMENTS 



In the following, embodiments of the 
present invention are described with reference to 
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the accompanying drawings. 

Fig. 3 is a block diagram showing a 
configuration example of a time-varying multi-path 
generating apparatus 9 according to an embodiment of 
5 the present invention. In the present specification, 
"time-varying" means a change and fluctuation of a 
characteristic along the time axis. For example, 
propagation path characteristics such as direction, 
phase, time delay, etc. are variable over time, i.e., 

10 are functions of time, with appropriate initial 
values assigned. The time-varying multi-path 
generating apparatus 9 according to the present 
embodiment includes an input device 1, a parameter 
input and control unit 2, a data storage unit 3, a 

15 random number generating unit 4, a propagation path 
generating unit 5, a time-varying function 
generating unit 6, an output data buffer memory unit 
7, and a time-varying propagation path output 
interface unit 8. Fig. 4 is a block diagram showing 

20 a configuration example of a multi-path simulator 14 
according to the embodiment of the present invention. 
The multi-path simulator 14 includes the time- 
varying propagation path output interface unit 8, 
the time-varying multi-path generating apparatus 9, 

25 a propagation path output unit 10, a signal 
synthesizing unit 11, a digital signal input 
terminal 12 , and a digital signal output terminal 13. 
Each of the items above is described below. 

The input device 1 is a device for 

30 inputting and outputting data files to/from the data 
storage unit 3, and inputting various parameters, 
required for the propagation path generation, that 
are displayed by the parameter input and control 
unit 2, and is constituted by a selection from a 

35 keyboard, a mouse, a floppy (registered trademark) 
disk, a memory card, a CD-ROM, a DVD, an MO disk, 
etc . 
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The parameter input and control unit 2 
displays information about propagation path 
generation parameters and various data files, and 
inputs required parameters using the input device 1. 
5 Here, data files containing information, such as a 
time-varying function (a function that varies over 
time) and an antenna directional-gain pattern, 
required for propagation path generation, are 
beforehand input through the parameter input and 

10 control unit 2, and stored in the data storage unit 
3 as required. Where the necessary data are 
beforehand stored in the data storage unit 3, it is 
not necessary to input the data through the 
parameter input and control unit 2. In this case, 

15 the parameter input and control unit 2 may only 
function as a parameter control unit. Further, 
parameters required for propagation path generation 
are specified, and the parameter input and control 
unit 2 controls all the propagation path generation 

20 steps, directing necessary information exchange 
between the data storage unit 3, the propagation 
path generating unit 5, and the time-varying 
function generating unit 6. 

In the data storage unit 3, a propagation 

25 path generation parameter file, a time-varying 

function file, an antenna directional-gain pattern 
file, an incoming propagation path direction 
density-function file, a distribution function file, 
etc . , are stored . 

30 The propagation path generation parameter 

file contains: 

a frequency of the propagation path F 

[GHz], 

the number N of propagation paths to be 

35 generated, 

the number M of output antenna branches, 
path parameter initial value generation 
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conditions , 

specification of a delay profile function, 
an average delay spread S D , 

a time-varying function type specification 

5 value; 

a constant-parameter generating condition 

file 

an average <T A > and standard deviation a A 
of shadow time interval T A [sec] , 

10 an average <L B > and standard deviation cr B , 

of shadow amplitude L B [dB] , 

an average <T C > and standard deviation ct c 
of shadow occurrence time interval T c .[sec] 

an average <C d > and standard deviation cr cd 
15 of falling coefficient C d , 

an average <C r > and standard deviation a Cr 
of rising coefficient C r , 

shadow time standard deviation a s (L B ) , 

inter -parameter correlation 
20 characteristics files p A (Q), p B (Q), and p c (^) 
containing inter-parameter correlation 
characteristics for an incoming propagation path 
angle difference, 

a shadow fluctuation calculation model 

25 file, 

a probability density function file 
containing probability density functions of shadow 
dimensions W h [m] , D h [m] , and H h [m] , 

a penetration angle probability density 
30 function file containing a probability density 

function of penetration angle 8 S [degrees], which is 
an angle of a shadowing object to the propagation 
path, 

a distance L 0 (the shortest propagation 
35 path) between the base station and the mobile 
station, 

a penetration position probability density 
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function file containing a probability density 
function of penetration position Ld [m] of the 
shadow, 

a moving speed probability density 
5 function file containing a probability density 

function of the moving speed Vs [m/s] of the shadow, 

moving speed V [m/s] of the mobile station, 
data sampling interval Dt [sec] , and the like, 

which are input from the parameter input 
10 and control unit 2. 

The time-varying function file contains 
plural shadow fluctuation functions Aj(t) and <j>j (t) 
that represent an amplitude and phase, respectively, 
of a propagation path specified by an input from the 
15 parameter input and control unit 2. 

The antenna directional-gain pattern file 
contains plural directional-gain patterns (V- 
polarization patterns and H-polarizat ion patterns) 
of various antennas to be used for simulation. The 
20 patterns may be obtained either by measurement or by 
calculation. If the patterns are not described in 
the same format, the patterns are converted into the 
same format through the parameter input and control 
unit 2 . 

25 The incoming propagation path direction 

density function file contains plural probability 
density functions that specify a propagation path 
occurrence probability by direction. One of the 
functions is selected according to environmental 

30 conditions when a propagation path is to be 
generated. 

The distribution function file contains 
probability density functions (e.g., an exponential 
distribution. function, an exponentiation 
35 distribution function, and a uniform distribution 
function) for generating random numbers 
corresponding to various environmental conditions. 
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The random number generating unit 4 
generates random numbers as required by the random 
number generation parameter that is specified by the 
parameter input and control unit 2, and provides the 
5 generated random numbers to the propagation path 
generating unit 5 and the time-varying function 
generating unit 6. When random number generation is 
required to continue by the parameter input and 
control unit 2, random numbers are repeatedly 

10 generated until a deactivation request is provided 
by the parameter input and control unit 2, and the 
generated random numbers are continuously provided 
to the propagation path generating unit 5 and the 
time-varying function generating unit 6. 

15 The propagation path generating unit 5 

carries out processes for generating propagation 
paths according to a generation request from the 
parameter input and control unit 2. The propagation 
path generating unit 5 first sets up initial values 

20 of the propagation path parameters (amplitude, phase, 
time delay, arrival direction, and the like) for N 
propagation paths based on the random numbers 
generated by the random number generating unit 4, 
and based on contents of the propagation path 

25 generation parameter file specified by the parameter 
input and control unit 2. Whether the propagation 
path generating unit 5 generates the time-varying 
shadow fluctuation based on the time-varying 
amplitude function, and the time-varying phase 

30 function stored in the data storage unit 3, or based 
on a time-varying amplitude function and a time- 
varying phase function that are to be calculated by 
the time-varying function generating unit 6, is 
determined according to contents of the propagation 

35 path generation parameter file. 

In the case that the time-varying function 
stored in the data storage unit 3 is to be used, the 
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time-varying function is read from the time-varying 
function file. Plural constant parameters for the 
time-varying function are generated based on the 
random numbers output by the random number 
5 generating unit 4 according to the contents of the 
propagation path generation parameter file, and N 
time-varying functions corresponding to N 
propagation paths are generated. 

In the case that the time-varying 

10 amplitude function and the time-varying phase 
function newly generated by the time-varying 
function generating unit 6 are to be used, the 
shadow fluctuation characteristics (amplitude and 
phase) that the time-varying function generating 

15 unit 6 calculates and outputs based on the contents 
the propagation path generation parameter file are 
used as the time-varying function, and N time- 
varying functions for the N propagation paths are 
generated with different shadow conditions that are 

20 obtained by using random numbers generated by the 
random number generating unit 4 based on the 
contents of the propagation path generation 
parameter file. 

The propagation path generating unit 5 

25 repeats the propagation path generating process, 

generating N time-varying functions every process, 
until a deactivation request is received. Every time 
a set of the N time-varying functions is generated, 
the set is cascaded such that a continuous time- 

30 varying function, including a continuous time- 
varying amplitude functions and a continuous time- 
varying phase function, for each of the N 
propagation paths is obtained. The continuous time- 
varying amplitude function is deducted from the 

35 initial amplitude value, and the continuous time- 
varying phase function is added to the initial phase 
value for each of the N propagation paths such that 
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t.he time-varying amplitude characteristics and the 
time-varying phase characteristics, respectively, 
are generated for each of the N propagation paths. 

As shown in Fig. 5, the propagation path 
5 generating unit 5 calculates phase rotation AQi(t) 
of each propagation path by the following Equation 1 
using information about the moving speed of the 
mobile station and about the direction of the 
initial incoming propagation path (specified in 

10 reference to the moving direction of the mobile 
station) , which phase rotation is added to the 
initial phase value of each propagation path such 
that the time-varying phase characteristics of the N 
propagation paths have a Doppler phase rotation 

15 effect. 

[Equation 1] 
AQi(t) = (V x t/X) • (cos(|>i x sinGi) 

20 Here, V [m/sec] represents the moving 

speed of the mobile station, t [sec] represents the 
time, Qi= (0i, <|>i) represents the initial arrival 
direction. of the i-th propagation path, and 
A, represents the wavelength of a wave coming through 

25 the propagation path. 

Furthermore, the propagation path 
generating unit 5 is capable of reading data of a 
specified antenna directional-gain pattern file from 
the data storage unit 3 through the parameter input 

30 and control unit 2, the file being selected based on 
the contents of the propagation path generation 
parameter file, generating propagation paths by 
multiplying the initial amplitude and the 
directional gain according to the initial arrival 

35 direction of each propagation path, and generating 
and outputting N time-varying propagation paths 
divided into M (M>=1) groups for each of <M antenna 
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patterns . 

The time-varying function generating unit 
6 calculates time-varying shadow characteristics 
(amplitude and phase) for the N propagation paths 
5 while varying calculation models and various 

calculation conditions for model calculations using 
random numbers generated by the random number 
generating unit 4 based on the contents of the 
propagation path generation parameter file specified 

10 by the parameter input and control unit 2, the 

calculation results of which being provided to the 
propagation path generating unit 5. The time-varying 
shadow functions are continuously calculated 
according to requests from the parameter input and 

15 control unit 2 or the propagation path generating 
unit 5 until a calculation deactivation request is 
provided. 

The output data buffer memory unit 7 
temporarily stores output data concerning the time- 

20 varying propagation paths generated by the 
propagation path generating unit 5. While a 
deactivation request is not provided by the 
parameter input and control unit 2, the stored data 
are provided to the output interface unit 8 when the 

25 volume of the data stored exceeds a predetermined 
threshold of data capacity. When a deactivation 
request from the parameter input and control unit 2 
is provided, all data that are stored, and data that 
are provided from the propagation path generating 

30 unit 5, are sequentially provided to the output 
interface unit 8 . 

With reference to Fig. 4, the output 
interface unit 8 serves as an interface unit for 
outputting the N time-varying propagation path 

35 characteristics (time-varying amplitudes, time- 
varying phases, initial time delays and initial 
arrival directions) generated by the propagation 
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path generating unit 5 to an external device such as 
the propagation path output unit 10. Data 
representing the time-varying propagation path 
characteristics are temporarily stored in the output 
5 data buffer memory unit 7, and are then provided to 
the output interface unit 8- The time-varying 
propagation path characteristics are output to the 
external device through the interface unit 8 as NxM 
time-serial values . 

10 The propagation path output unit 10 

receives the N time-varying propagation paths 
calculated and generated by the propagation path 
generating unit 5 from the output interface unit 8, 
and converts the N time-varying propagation paths 

15 into real parts and imaginary parts of a complex 
amplitude using N time-varying amplitudes and N 
time-varying phases; and then, carries out digital- 
to-analog conversion to obtain and output 2xNxM 
analog signals, 

20 The signal synthesizing unit 11: 

receives M groups of time-varying 
propagation path signals, the groups consisting of 
the real part time-varying signals W R {11) through 
Jtf R (MN) of the N complex amplitudes, and the imaginary 

25 part time-variable signals W T {11) through Jtfj <MN) of the 
N complex amplitudes, 

distributes an IF or RF signal u (t) to M 
antenna branches, the IF or RF signal being input to 
the digital input signal terminal 12, 

30 generates orthogonal signals (Q component) , 

the phases of which are shifted by 7i/2 from the 
distributed signals for every digital signal, 

each of the digital input signals, which 
is the I component u T (t), and the Q component u Q (t) 
35 is input into transversal type circuits (refer to 
Fig. 16) , each of which consists of cascade 
connections of N-l delay elements, 
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N signals of the I component having 
different time delays corresponding to time delay of 
each of the propagation paths provided by the 
propagation path output unit 10 and corresponding 
5 real part time-varying signal of the N complex 
amplitudes are multiplied, and similarly, the N 
signals of the Q component having different time 
delays corresponding to time delay of each of the 
propagation paths provided by the propagation path 

10 output unit 10 and corresponding imaginary part 

time-varying signal of the N complex amplitudes are 
multiplied. The multiplied results are added for 
each of the N propagation paths. This process is 
repeated M times such that M signals y (M> (t) are 

15 generated, and are output from the digital signal 
output terminal 13. Equation 2 that follows shows 
the relations between the input signal and the 
output signal. 

[Equation 2] 

y(M)(t) = yi(M)(t) + y Q <M)(t) 

20 

In Equation 2, the carrier component is 
omitted from notation. 

Fig. 6 is a flow chart showing operations 
of the time-varying multi-path generating apparatus 

25 9 according to the embodiment of the present 

invention. The operations of the present invention 
are explained in reference to steps SI through S15 
shown in the flow chart. Here, it is presumed that 
the data files, such as the propagation path 

30 generation parameter file, the time-varying function 
file, the antenna' directional-gain pattern file, and 
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the incoming propagation path direction density 
function file are stored in the data storage unit 3 
by an operator in advance of generating time-varying 
multi-path propagation paths according to the 
5 present invention. 

At step SI, the parameter input and 
control unit 2 reads the propagation path generation 
parameter file from the data storage unit 3, and 
secures memory capacity in the output data buffer 

10 memory unit 7, the capacity being sufficient for 
storing data that are to be generated for N 
propagation paths and M output antenna branches. The 
parameter input and control unit 2 shares the 
propagation path generation parameter file with the 

15 propagation path generating unit 5, or provides the 
file to the propagation path generating unit 5, and 
proceeds to step S2 . 

At step S2, the parameter input and 
control unit 2 reads M antenna directional-gain 

20 pattern files and direction density function files 
of an incoming propagation path corresponding to M 
output antenna branches from the data storage unit 3, 
and proceeds to step S3. 

At step S3, the propagation path 

25 generating unit 5 requests the random number 

generating unit 4 to provide random numbers that 
fulfill the propagation path parameter initial value 
generation conditions of a propagation path to be 
generated, and obtains the random numbers from the 

30 random number generating unit 4. The random numbers 
are used to set up the initial conditions (initial 
amplitude, initial phase, initial time delay, and 
initial arrival direction) of N propagation path 
parameters. 

35 As for the initial amplitude and the 

initial time delay, the propagation path generation 
parameter file may specify an exponential 
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distribution function with a predetermined 
coefficient to be the delay profile function. In 
this case, the initial amplitude and initial time 
delay of N propagation paths are determined by a 
5 . delay time step AT, the number N of propagation 

paths to be generated, and an average delay spread 
value S D of a propagation path to be generated. • 
Further, if the initial phase is uniformly 
distributed in the range between 0 and 2n, the 

10 initial phase is determined by generating N uniform 
random numbers; and the initial arrival direction is 
determined by generating N random numbers according 
to the direction density function specified by the 
direction density function file of the incoming 

15 propagation path. How the random numbers for 

generating the initial values are to be generated is 
prescribed by the propagation path parameter initial 
value generation conditions specified by the 
propagation path generation parameter file. 

20 After the initial values for the N 

propagation path parameters are determined, the 
initial amplitude values and the antenna directional 
gain in the direction of each incoming propagation 
path read from the antenna directional-gain pattern 

25 file stored in the data storage unit 3 are 

multiplied such that NxM propagation path parameters 
for N propagation paths and M antenna branches are 
set up. Then, the process proceeds to step S4. 

Fig. 7 shows examples of initial values of 

30 the propagation path parameters generated according 
to the propagation path parameter initial value 
generating conditions, wherein the number of the 
propagation paths is set at N, the shape of the 
delay profile is set rectangular (i.e., amplitudes 

35 are constant regardless of delay time) , the delay 

step is set at 1 jisec, the arrival direction is set 
at a uniform distribution over a 0 - 2n range, the 



-23- 



phase characteristics are set at a uniform 
distribution over the 0 - 2n range, and the moving 
speed of the mobile station is set at 0 m/sec. At 
section (a) of Fig. 7, the delay profile constituted 
5 by the initial time delay and initial amplitude is 
shown. At section (b) of Fig. 7, the initial phase 
characteristics are shown. At section (c) of Fig. 7, 
the initial arrival direction characteristics are 
shown. The propagation path parameter initial value 

10 generation conditions of a propagation path to be 

generated are determined in the manner as described 
above, for example. 

At step S4, the parameter input and 
control unit 2 determines whether a function to be 

15 generated by model calculation carried out by the 
time-varying space generating unit 6 is to be used 
as the time-varying function based on a type 
specification value for the time-varying function 
contained in the propagation path generation 

20 parameter file read from the data storage unit 3. 

At step S5 (step S4 NO) , the parameter 
input and control unit 2 reads the time-varying 
function file from the data storage unit 3, and 
shares the file with the propagation path generating 

25 unit 5, or provides the file to the propagation path 
generating unit 5. Based on the time-varying 
function file, the type of the time-varying function 
is determined, and the constant parameter generation 
condition file of the time-varying function is read 

30 from the propagation path generation parameter file 
stored in the data storage unit 3. Further, the data 
sampling interval is read, and the inter-parameter 
correlation characteristic file concerning a 
propagation path angle difference is also read, and 

35 the process proceeds to step S6. 

At step S6, the propagation path 
generating unit 5 requires the random number 
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generating unit 4 to generate and output random 
numbers corresponding to the probability density 
functions of the shadow time interval T A/ the shadow 
amplitude L B , and the shadow occurrence time 
5 interval T c contained in the propagation path 

generation parameter file, and obtains the random 
numbers. If correlations are required to be attached 
to these random numbers, correlations corresponding 
to the difference of an incoming propagation path 

10 angle from the initial arrival direction specified 

by the parameter correlation characteristic file are 
attached such that random numbers having the 
requested correlations are generated. In this manner, 
the constant parameters of the time-varying function 

15 are determined using the random numbers, and the 
process proceeds to step S7. 

Here, an example of a time-varying 
amplitude function serving as a time-varying 
function file is explained with reference to graphs 

20 presented in Fig. 8, which graphs show a shadow 

effect of a human being and the like. As shown at 
section (a) of Fig. 8, the shadow effect is 
represented by a shadow amplitude L B [dB] that is a 
function of time f(t), wherein T A [sec] represents a 

25 shadow time interval, and T c [sec] represents a 
shadow occurrence time interval. Since T A and T c 
determine the shape of the function as shown at 
section (b) of Fig. 8, T A and T c serve as constant 
parameters. Then, as conditions that determine these 

30 constant parameters, the following is premised, for 
example. Namely, the shadow time interval T A [sec] 
follows a normal distribution having an average <T A > 
[sec] and a standard deviation a A [sec] ; the shadow 
amplitude L B [dB] follows a normal distribution 

35 having an average <L B > [dB] and a standard deviation 
[dB] ; and the shadow occurrence time interval T c 
[sec] follows a normal distribution having an 



average <T C > [sec] and standard deviation a c [sec] - 
Then, each of the parameters T A and T c can be 
obtained by N random numbers according to 
distribution of each of the parameters, such random 
numbers being required, and generated by the random 
number generating unit 4. 

Fig. 9 shows another example of the time- 
varying amplitude function serving as the time- 
varying "function file, when a big obstacle, such as 
a truck, or a bus causes the shadow effect. At 
section (a) of Fig. 9, the shadow effect is 
expressed as a shadow amplitude . L B [dB] that is a 
function of time g(t), where T A [sec] represents a 
shadow time interval, T c [sec] represents a shadow 
occurrence time interval, C d [dB/sec] represents a 
falling coefficient, C r [dB/sec] represents a rising 
coefficient, and a s [dB] represents a standard 
deviation of shadow amplitude L B [dB] concerning the 
time of shadow occurrence. Since the function g(t) 
takes various shapes depending on these values as 
shown at section (b) of Fig. 9, these values are 
constant parameters that determine the shape of the 
function. Then, concerning conditions that determine 
these constant parameters, it is premised as follows 
for example. Namely, the shadow time interval T A 

[sec] follows a normal distribution having an 
average <T A > [sec] and a standard deviation a A [sec] 
The shadow amplitude L B [dB] follows a normal 
distribution having an average <L B > [dB] and a 
standard deviation a B [dB] . The shadow occurrence 
time interval T c [sec] follows a normal distribution 
having an average <T C > [sec] and a standard 
deviation a c [sec] . The falling coefficient C d 

[dB/sec] follows a normal distribution having an 
average <C d > [dB/sec] and a standard deviation a cd 

[dB/sec] . The rising coefficient C r [dB/sec] follows 
a normal distribution having an average <C r > 
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[dB/sec] and a standard deviation a Cr [dB/sec] . The 
shadow fluctuation follows a normal distribution 
. having an average 0 and a standard deviation a s [dB] . 
Then, normally distributed random numbers serve each 
5 distribution, the random number generating unit 4 
generating N random numbers accordingly. The time- 
variable functions that represent cases such as 
shown in Fig. 8 (shadow by a human) and Fig. 9 
(shadow by a big obstacle such as a truck or a bus) 

10 according to predetermined environmental conditions 
as required are beforehand prepared, and stored in 
the data storage unit 3 such that desired 
propagation path simulations are performed. 

At step S7, based on the constant 

15 parameters of the time-varying functions determined 
at Step 6, the propagation path generating unit 5 
generates time-varying shadow amplitude 
characteristics and time-varying shadow phase 
characteristics of N propagation paths, the time- 

20 varying shadow amplitude characteristics are 
deducted from the initial amplitude of each 
propagation path, and the time-varying shadow phase 
characteristics are added to the initial phase of 
each propagation path such that shadow fluctuation 

25 characteristics are generated for the amplitude and 
the phase, respectively, of each propagation path. 
At the same time, the propagation path generating 
unit 5 calculates the time phase rotation AQi(t) by 
the formula (1) using information about the moving 

30 speed of the mobile station, and information about 
the incoming direction of an initial propagation 
path (specified in reference to the moving direction 
of the mobile station) , which time phase rotation is 
added to the initial phase of each propagation path 

35 such that a Doppler phase rotation effect can be 

attached to the time-varying phase characteristics 
of each of the N propagation paths. Then, the 
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process proceeds to step S8. 

At step S8, the propagation path 
generating unit 5 confirms that the time-varying 
data that are obtained by the next random number 
5 generation and that are to be added in continuation 
would not influence the present result, then the 
data are transferred to the output data buffer 
memory 7. Then, the process proceeds to step S9. 

At step S9, while the propagation path 

10 generating unit 5 repeatedly generates data of the 

shadow fluctuation characteristics for the amplitude 
and the phase of each propagation path, and 
repeatedly adds serially in the time domain such 
that output data are generated, the propagation path 

15 generating unit 5 determines whether there is a 
propagation path generation deactivation request 
from the parameter input and control unit 2. If it 
is determined that no deactivation requests are 
present, the process returns to step S6, and 

20 generation of the shadow fluctuation characteristics 
for the amplitude and the phase of subsequent 
propagation paths is repeated. Here, the time- 
varying shadow amplitude characteristics to be 
deducted and the time-varying shadow phase 

25 characteristics to be added for each propagation 
path are continuously generated in every shadow 
occurrence time interval. 

Therefore, when the shadow time interval 
T A [sec] is greater than the shadow occurrence time 

30 interval T c [sec], shadows overlap as shown by the 
second graph in Fig. 10. In this case, the shadow 
amplitudes of the shadows are simply added such that 
the time-varying shadow amplitude characteristics 
and the time-varying shadow phase characteristic are 

35 obtained as shown by the bottom graph in Fig. 10. 

In the case that the parameter input and 
control unit 2 provides a propagation path 
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generation deactivation request, the process proceed 
to step S 9 A wherein data generation is stopped, data 
that have been generated before the deactivation 
request are provided to the output data buffer 
5 memory 7, and the process proceeds to step S10. 

At step S10, the output data buffer memory 
unit 7 stores the data provided by the propagation 
path generating unit 5, and outputs the data to the 
output interface unit 8 when the amount of the data 

10 exceeds a predetermined value. In the case that a 

generation deactivation request is provided from the 
parameter input and control unit 2, all the data 
stored are sequentially output to the output 
interface unit 8 until all the data are output, even 

15 if the data amount of the output data buffer memory 
unit 7 does not reach the predetermined value. In 
this manner, the propagation path generation is 
finished. 

At step Sll (step S4 YES), the parameter 
20 input and control unit 2 determines which 

calculation model is to be used by the time-varying 
function generating unit 4 based on the shadow 
fluctuation calculation-model file, frequency F 
[GHz], and the- distance L 0 [m] between the base 
25- station and the mobile station (the shortest of the 
propagation paths) . Then, the process proceeds to 
step S12. 

At step S12, the propagation path 
generating unit 5 requests the random number 

30 generating unit 4 to generate random numbers for 

each of dimensions (W h [m] , D h [m] , and H h [m] ) of a 
shadowing object (obstacle) , the penetration angle 
0s [degrees] of the shadowing object to the 
propagation path, the penetration position Ld [m] of 

35 the shadowing object, and the moving speed Vs [m/s] 
of the shadowing object based on the probability 
density function of the shadowing object dimensions 
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(W h [m] , D h [m] , H h [m] ) , the penetration angle 0s 
[degrees] of the shadowing object to the propagation 
path, the probability density function of the 
penetration position Ld [m] of the shadowing object, 
5 and the probability density function of the moving 

speed Vs [m/s] of the shadowing object, respectively. 
The random numbers obtained in this manner are 
provided to the time-varying function generating 
unit 6, various parameters for the calculation model 

10 are determined, and the time-varying function 

generating unit 6 generates N time-varying functions. 
When correlations are required between the random 
numbers, a correlation corresponding to the 
difference between the initial arrival direction 

15 specified by the parameter correlation 

characteristic file and the incoming propagation 
path angle is attached, and the random numbers are 
generated. Then, the process proceeds to step S13. 

Fig. 11 is a flowchart of a calculation 

20 model for generating a time-varying amplitude 
function and a time-varying phase function. 

The outline of the calculation model is 
illustrated in Fig. 12. In the calculation model of 
the present embodiment shown in Fig. 12, Lo [m] 

25 represents the distance between the base station and 
the mobile station (the shortest propagation path) , 
F [GHz] represents frequency of the wave coming 
through the propagation path, Ld [m] represents the 
penetration position of the shadowing object 

30 (obstacle) measured from the mobile station, 

0 [degrees] represents the penetration angle of the 
shadowing object, and the dimensions of the 
shadowing object are given as an elliptic cylinder 
having a minor axis a [m] and a major axis b [m] . 

35 When the shadowing object is present right in the 

line of sight, the received electric field strength 
E is expressed as follows. 
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[Equation 3] 



E = E l +E 2 



(2) 



Here, El and E2 are the field strengths of 
the radio propagation paths that are diffracted by 
5 knife-edges 1 and 2, respectively. 

Further, El and E2 are given by the 
following formula . 

[Equation 4] 



strength at the mobile station, i.e., the receiving 
point where no shadows are present. Ef is expressed 
as follows . 



5S$x = l or 2 



(3) 



10 



Here, E f is the free-space electric field 



[Equation 5] 




-J30P 



[V/m] 



15 



Here, P represents the radiation power at 
a transmitting antenna. C and S are Fresnel 
integrals and are expressed as. follows. 

[Equation 6] 
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CifO x ) - jS(fO x ) = e- Ji *' 2)i2 dt >&x = 1 or 2 ( 4 ) 



Here, co is called a clearance function, 
and is expressed as follows. 
[Equation 7] 



Here, 1 0 [m] , equivalent to L 0 elsewhere, 
represents the distance between the base station and 
the mobile station, X [m] represents the wavelength 
of the wave coming through the propagation path, and 
10 (xl, yl) and (x2, y2 ) represent the coordinates of 

the knife-edges 1 and 2, respectively, shown in Fig. 
12. 

With reference to Fig. 11, the calculation 
process of the time-varying function based on this 

15 calculation model is explained. At step S401, the 
positions (coordinates) of the knife-edges 1 and 2 
in the case that the center of the shadowing object 
(obstacle) is located on the line-of -sight are 
calculated using the parameters. Namely, the 

20 parameters are the distance L 0 [m] between the base 

station and the mobile station (shortest propagation 
path) , ^nd the penetration position Ld [m] of the 
shadowing object measured from the mobile station, 
the shadowing object penetration angle 0 [degrees], 

25 and the shadowing object dimensions (in the present 
model, the elliptic cylinder having the minor axis a 
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[m] and the major axis b [m] ) . Then, the electric 
field strengths El and E2 of the propagation paths 
diffracted by the knife-edges 1 and 2, respectively, 
are calculated, and the electric field strength E is 
5 obtained by calculating a vector sum of El and E2 . 

At step S402, the position of the 
shadowing object is moved (downward in the case of 
Fig. 12) away from the line-of -sight by a step 
distance that is a product of the moving speed [m] 

10 of the shadowing object and a data sampling interval, 
the positions of the knife-edges 1 and 2 are 
determined like in step S401, the electric field 
strengths of the diffraction propagation paths 1 and 
2 are calculated, and the electric field strength E 

15 is obtained by calculating a vector sum of El and E2 . 
This process is r.epeated until the vector sum E 
exceeds E f that is the electric field strength where 
no shadows are present (electric field strength of 
the direct propagation path) . 

20 At step S403, the position of the 

shadowing object is moved in the direction opposite 
to the case described above (i.e., upward in the 
case of Fig. 12) away from the line-of -sight by the 
step distance that is the product of the moving 

25 speed [m] of the shadowing object and the data 

sampling interval, the positions of the knife-edges 
1 and 2 are determined like in step S401, the 
electric field strengths of the diffraction 
propagation paths 1 and 2 are calculated, and the 

30 electric field strength E is obtained by calculating 
the vector sum of El and E2 . This process is 
repeated until the vector sum E exceeds E f that is 
the electric field strength where no shadows are 
present (electric field strength of the direct 

35 propagation path) . 

At step S404, the amplitudes and the 
phases of the electric field strength E obtained at 
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steps S401, S402 and S403 are stored in the memory 
storage as the time-varying functions corresponding 
to the moving time of the shadowing object. 

In Fig. 13, a solid line curve illustrates 
a calculation example of the received-power 
fluctuation with the parameter value being set as 
follows: L 0 =13 . 5 m, F=7 GHz, Ld=3 m, 9=45 deg, 
a=0.49 m, and b=0.20 m. Here, the relations between 
the received power P and the electric field E are 
defined by the following formula. 

[Equation 8] 

P = 10 logio I E/Er 1 2 ( 

A dotted line curve in Fig. 13 represents 
a measurement result performed on the same 
conditions. In the measurement, a human body serves 
as the shadowing object. The horizontal axis 
represents the position of the human body in 
reference to the line-of -sight . As seen from Fig. 13 
the calculation model well approximates the 
received-power fluctuation, i.e., the time-varying 
amplitude function. 

With reference back to Fig. 6, at step S13 
the propagation path generating unit 5 time-serially 
(serially in the time domain) generates the shadow 
fluctuation characteristics for the amplitude and 
the phase of each propagation path by subtracting 
the time-varying shadow amplitude characteristics 
from the corresponding initial value, and by adding 
the time-varying shadow phase characteristics to the 
corresponding initial value, given that N time- 
varying functions (the time-varying amplitude 
functions and the time-varying phase functions) 
generated by the time-varying function generating 
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unit 6 represent the time-varying shadow amplitude 
characteristics and the time-varying shadow phase 
characteristics, respectively, of each of the N 
propagation paths . 
5 At this time, the propagation path 

generating unit 5 calculates the time phase rotation 
AQi(t) of each propagation path according to the 
formula (1) based on the moving speed of the mobile 
station and the information about the initial 

10 direction of the incoming propagation path 

(specified in reference to the moving direction of 
the mobile station) , and the time phase rotation is 
added to the initial phase of each propagation path 
such that the Doppler effect is attached to the 

15 time-varying phase characteristics of the N 

propagation paths. Then, the process proceeds to 
step S14. 

At step S14, the propagation path 
generating unit 5 starts outputting the. data 

20 generated at step S13 to the output data buffer 

memory 7 unless it is determined that time-varying 
data that are to be continuously added by the next 
random number generation influence the present 
result, while the shadow fluctuation characteristic 

25 data of the amplitude and the phase for each 

propagation path continuously generated at step 13 
based on the random numbers generated at step S12 
are repeatedly added in the time domain at step S15 
as shown in Fig. 10. Then, . the process proceeds to 

30 step S15 . 

At step S15, the propagation path 
generating unit 5 determines whether there is a 
propagation path generation deactivation request 
from the parameter input and control unit 2. If 
35 there is no deactivation request (step S15 NO) , the 
process returns to step S12, and the propagation 
path generating unit 5 repeats generation of the 
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shadow fluctuation characteristics for the amplitude 
and the phase of each propagation path. Here, the 
time-varying shadow amplitude characteristics and 
the time-varying shadow phase characteristics that 
5 are deducted and added, respectively, from/to the 
amplitude and the phase, respectively, of each 
propagation path are continuously generated for 
every shadow occurrence time interval. For this 
reason, a group of time-varying characteristics that 

10 overlaps with another group of time-varying 

characteristic may be generated as shown in Fig. 10 
as described in reference to step S9. If such 
overlap takes place, the overlapping shadow 
characteristics are simply added to obtain the time- 

15 varying shadow phase characteristics and the shadow 
amplitude characteristics. 

In the case where there is a propagation 
path generation deactivation request from the 
parameter input and control unit 2 (step S15 YES), 

20 the process proceeds to step S15A wherein data 
having been so far generated are sequentially 
provided to the output data buffer memory 7, and the 
process proceeds to step S10. 

Fig. 14 shows an example of propagation 

25 path status of 10 propagation paths after level 
fluctuations are attached. The example shows the 
time-varying amplitude fluctuation of each 
propagation path corresponding to the same 
penetration position, penetration angle, dimensions 

30 of the shadow, and moving speed of the shadow, where 
the moving speed of the mobile station is set at 0 
m/s, and the shadow occurrence interval is set by an 
exponential distribution having an average of 0.5 
sec. In the example graphs, the vertical axis 

35 represents the shadow amplitude L B [dB] , and the 

horizontal axis represents the time where the shadow 
time interval T A of each shadow is set at 1 sec. 
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Further, the example is prepared on the premise that 
no shadows overlap. Note that shadows overlap in the 
case of Fig. 10. Nevertheless, the conditions giving 
the shadow fluctuations are not limited to those 
5 described above, but are replaced with, or added to 
the above calculation conditions according to 
environmental conditions under which a simulation is 
to be performed. 

Fig. 15 shows a configuration example of 
10 the propagation path output unit 10 of the multi- 
path simulator 14 shown in Fig. 4 according to the 
embodiment of the present invention. The propagation 
path output unit 10 includes a real part /imaginary 
part separation and conversion unit 15, and digital- 
is to-analog (D/A) converters 16. The real 

part/imaginary part separation and conversion unit 
15 extracts a real part signal and an imaginary part 
signal from each of the MxN time-varying complex 
amplitudes of the propagation path fluctuations 
20 provided by the multi-path generating apparatus 9, 
and outputs the real part signals and the imaginary 
part signals to the D/A converters 16. There are 
2xMxN D/A converters 16 that carry out D/A 
conversion and output the time-varying propagation 
25 path signals in analog form. 

Fig. 16 shows a configuration example of 
the signal synthesizing unit 11 of the multi-path 
simulator 14 according to the embodiment of the 
present invention. The signal synthesizing unit 11 
30 includes an orthogonal signal generating unit 17, 
transversal type circuits 18, delay elements 19, 
multipliers 20, and adders 21. The analog output of 
the propagation path output unit 10 is provided to 
each of the multipliers 20 as a multiplication 
35 weight. The orthogonal signal generating unit 17 
generates an orthogonal signal (Q component) from 
the signal (I component) provided through the 
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digital signal input terminal 12. The I component 
and the Q component are divided into M signals that 
are named the first signal through the Mth signal. 
Here, N represents the number of propagation paths 
5 and M represents the number of antennas. The first 
signal of the I component is provided to the 
transversal type circuit 18 numbered as 1-1, the 
second signal is provided to the transversal type 
circuit 18 numbered as 1-2, and so on up to the Mth 

10 signal being provided to the transversal type 

circuit 18 numbered as I-M. Similarly, the first 
signal of the Q component is provided to the 
transversal type circuit 18 numbered as Q-l, and the 
second signal is provided to the transversal type 

15 circuit 18 numbered as Q-2, and so on up to the Mth 
signal being provided to the transversal type 
circuit 18 numbered as Q-M. 

As described above, the propagation path 
output unit 10 outputs N propagation path signals 

20 for each of M antennas, the signals having a real 
part signal and an imaginary part signal. Namely, 
each of the real part signals are called wR (11> , 
wR (12) , and through wR (1N) ; wR (21) , wR (22) , and through 
wR (2N) ; and up through wR (M1) , wR (M2) , and through wR (MN) . 

25 Similarly, the imaginary part signals are called 
wl (11) , wl (12) , and through wI (1N) ; wl (21) , wl (22) , and 
through wI (2N) ; and up through wI (M1) , wI (M2) , and 
through wI (MN) . 

The propagation path fluctuation signals 

30 wR (11> , wR (12> , and through wR (1N) serve as the weights 
provided to the first, second, and through Nth 
multipliers 20, respectively, of the transversal 
type circuit 18' No. 1-1. Similarly, the propagation 
path fluctuation signals wR (M1) , wR (M2) , and through 

35 wR (MN) serve as the weights provided to the first, 

second and through Nth multipliers 20, respectively, 
of the transversal type circuit 18 No. I-M . 
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Similarly, the propagation path 
fluctuation signals wl (11) , wl (12) / and through wI (1N) 
serve as the weights provided to the first, second, 
and through Nth multipliers 20, respectively, of the 
5 transversal type circuit 18 No. Q-l; and the 

propagation path fluctuation signals wI (M1) , wI (M2) , 
and through wI (MN) serve as the weights provided to 
the first, second, and through Nth multipliers 20, 
respectively, of the transversal type circuit 18 No. 

10 Q-M. The I component and Q component output by each 
transversal type circuit 18 are added for every 
group corresponding to each of the M antennas, and M 
output signals are obtained.. In this manner, the 
input signal is transformed into the output signals 

15 wherein the influence of a directional pattern of 
every antenna is taken into consideration. 

In Fig. 17, thin lines show examples of 
outputs of the fading simulator according to the 
embodiment of the present invention. Thick lines 

20 show Rayleigh distributions where a parameter K 

takes values such as K=-infinity, 3, 6, 10, and 15 
dB, respectively. 

Further, in Fig. 18, thin lines show 
examples of measurement results of the receiving 

25 level characteristics when measured inside of a 

building. Thick lines show Rayleigh distributions 
where K=-infinity, 6, and 10 dB, respectively. As 
seen from Fig. 17 and Fig. 18, both examples of the 
output results of the fading simulator according to 

30 the embodiment and the measurement results are well 
in agreement with the Rayleigh distributions with 
various K factors. 

[Effect of the Invention] 
As described above, according to the 

35 present invention, the time-varying characteristics 
(the time-varying multi-path fluctuation) of the 
propagation path parameters (power level 
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characteristics, phase characteristics, time delay 
characteristics, and arrival direction 
characteristics) corresponding to desired simulation 
conditions are generated by synthesizing propagation 
5 path parameters of individual propagation paths, 
taking correlations of spatial propagation paths 
into consideration, rather than simply employing 
conventional Rayleigh fading or Nakagami-Rice fading. 
Further, it was confirmed that Rayleigh fading, 

10 Nakagami-Rice fading, and the like were reproducible, 
even when the mobile station is not moving, by 
obtaining a vector sum of two or more propagation 
paths generated according to the embodiment of the 
present invention. In this manner, the multi-path 

15 fading simulator that can simulate propagation path 
fluctuations that occur by actual propagation 
mechanisms is realized, which is highly valid and 
important . 

Further, the present invention is not 
20 limited to these embodiments, but various variations 
and modifications may be made without departing from 
the scope of the present invention. 

The present application is based on 
Japanese Priority Application No. 2003-045104 filed 
25 on February 21, 2003 with the Japanese Patent Office, 
the entire contents of that are hereby incorporated 
by reference. 



